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Abstract. The interaction of an expanding Pair- 
Electromagnetic pulse (PEM pulse) with a shell of bary- 
onic matter surrounding a Black Hole with electromag- 
netic structure (EMBH) is analyzed for selected values 
of the baryonic mass at selected distances well outside 
the dyadosphere of an EMBH. The dyadosphere, the re- 
gion in which a super critical field exists for the creation 
of e+e^ pairs, is here considered in the special case of 
a Reissner-Nordstrom geometry. The interaction of the 
PEM pulse with the baryonic matter is described using a 
simplified model of a slab of constant thickness in the lab- 
oratory frame (constant-thickness approximation) as well 
as performing the integration of the general relativistic hy- 
drodynamical equations. The validation of the constant- 
thickness approximation, already presented in a previous 
paper Ruffini, et al.(1999) for a PEM pulse in vacuum, is 
here generalized to the presence of baryonic matter. It is 
found that for a baryonic shell of mass-energy less than 
1% of the total energy of the dyadosphere, the constant- 
thickness approximation is in excellent agreement with 
full general relativistic computations. The approximation 
breaks down for larger values of the baryonic shell mass, 
however such cases are of less interest for observed Gamma 
Ray Bursts (GRBs). On the basis of numerical computa- 
tions of the slab model for PEM pulses, we describe (i) the 
properties of relativistic evolution of a PEM pulse collid- 
ing with a baryonic shell; (ii) the details of the expected 
emission energy and observed temperature of the associ- 
ated GRBs for a given value of the EMBH mass; IO'^Mq, 
and for baryonic mass-energies in the range 10^'* to 10^^ 
the total energy of the dyadosphere. 

Key words: black holes physics - gamma rays: bursts - 
gamma rays: theory - gamma rays: observations 



1. Introduction 

That vacuum polarization process a la Heisenberg-Euler- 
Schwinger ( Heisenberg fc Euler 1931 , Bchwinger 1951 ) 



can occur in the field of a Kerr Newmann EMBH and 
that they naturally lead to a model for gamma-ray bursts 



was pointed out in Damour & Ruffini (1975). The basic 



energetics of the process, governed by the Christodoulou- 
Rufhni mass formula, for an E MBH gives as shown in 
IChristodoulou fc Ruffini (1971) 



M 



) = 167r 



with 

-l(5)(g^ + 4.V)<i, 



(1) 

(2) 
(3) 



where M, L and Q are respectively the mass energy, the 
angular momentum and the charge of the EMBH and Mjr 
is the irreducible mass, r+ is the horizon radius, is the 
quasi-spheroidal cylindrical coordinate of the horizon eval- 
uated at the equatorial plane and S is the horizon surface 
area. Extreme black holes satisfy the equality in Eq.(||). 
The vacuum polarization process being "reversible" trans- 
formations in the sense of IChristodoulou & Ruffini (1971) 



can extract an energy up to 29% of the mass-energy of an 
extremal rotating black hole and up to 50% of the mass- 
energy of an extremely magnetized and charged black hole. 

Although in general such a process is endowed with 
axial symmetry, in order to clarify the pure interplay of 
the gravitational and electrodynamical phenomena and 
also for simplicity, we have examined (Ruffini 1998 and 



Preparata et al. 1998a, t) the limiting cases of such a pro- 



cess in the field of a Reissner-Nordstrom geometry whose 
spherically symmetric metric is given by 

(fs = 9tt{r)cPt + grr{r)d?r + r'^SO + sin^ ed?(t) , (4) 



where gtt{r) 



1 



2GM 



-a'^{r) and 



grrjr) = a ^(r). The dyadosphere, defined (see Fig. 
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1 of Preparata et al. 1998a, b) as the region where the 

2 3 

electric field exceeds the critical value, £r — "I 
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(Heisenberg & Euler 1931, Schwinger 1951), where m and 
e are the mass and charge of the electron, extends between 
the horizon radius 



= 1.47 • 10'V(1 + Vl-C^) cm. 



(5) 



where we have introduced the dimcnsionless mass and 



charge parameters /i — 
outer radius 

T-ds = 1.12 • 10*VM?cm. 



M 



(mVg) 



< 1, and the 



(6) 



Consequently the local number density of electron and 
positron pairs created in the dyadosphere as a function of 
radius 



Q 



1 



47rr2 (^) e 

V mc / 

and their energy density is given by 
n2 / / v4 



87rr4 



1 



(7) 



(8) 



(see Figs. [2] & [3] of [Prcparata et al. 1998a,b|) . The total 
energy of pairs converted from the static electric energy 
and deposited within the dyadosphere is then 



1- ( — 

rds 



and the total number of electron and position pairs in the 
dyadosphere is 



N°+ 

e+e 



Q-Qc 



1 



[rds 



-) 



h 

mc 



(10) 



where Qc = £cf'\ (sec [Prcparata et al. 1998a,b ). In the 
hmit of ^ 0, Eq.(|) leads to Edy^, which 
coincides with the energy extractable from EMBHs by 
reversible processes (Mir = const.), namely E — il/jr = 
i^(see Fig. 4 of [Preparata et al. 1998a,b[ ). 

Due to the very large pair density given by Eq. (^ and 
to the sizes of the cross-sections for the process e+e^ ^ 
7 -|- 7, the system is expected to thermalize to a plasma 
configuration for which 



(11) 



where N\ _ is the number of e^e -pairs created in the 



dyadosphere (see Preparata et al. 1998a, b). These are the 
initial conditions for the evolution of the dyadosphere. In a 



previous paper (Ruffini et al. 1999) we presented the tem- 
poral evolution of the dyadosphere in vacuum giving origin 
to an extremely sharply defined and extremely relativistic 
expanding pulse of pair and electromagnetical radiation of 
a constant length in the laboratory frame: the PEM pulse. 
In this paper we present results of the collision of the PEM 
pulse with a remnant of baryonic matter surrounding the 



just formed black hole. We assume the PEM pulse to col- 
lide with a shell of baryonic matter of constant density 
and at a radius of the order of 100 times the radius of the 
dyadosphere rds Eq.(|^). The shells have this thickness of 
the order of 10 times rds- The mass-energies of the shells 
are taken to be 10^^ — 10^^ of the total energy of the 
dyadosphere. 

The work contains the following sections: Sect. || 
presents discussions of the hydrodynamical equations of 
a PEM pulse interacting with the baryonic shell; Sect. 
^ defines the parameters of the baryonic shells, the be- 
haviour of the PEM pulse colliding with a baryon shell 
as well as before and after the collision are presented; in 
Sect. ^ the results of a numerical calculation solving the 
hydrodynamical equations of Sect. § are compared to the 
results of the analytical model of slab approximation for 
selected parameters. 

2. General relativistic hydrodynamical equations for 
PEM pulse in presence of baryonic matter 

The evolution of the PEM pulse in vacuum was treated 
in a previous work (Ruffini et al. 1999). We here gener- 
alize that treatment to the case in which baryonic mat- 
ter is present and we outline the relevant relativistic 
hydrodynamic equations. As in the previous treatment 



(9) (Ruffini et al. 199£), we assume the plasma fluid of 



pairs, photons and baryonic matter to be described by the 
stress-energy tensor 



Tf"" = pg'"' + {p + p)U''U'' 



(12) 



where p and p are respectively the total proper energy 
density and pressure in the comoving frame of the plasma 
fluid. The [/^ is the four-velocity of the plasma fluid. We 
have 



r\2 



gttiUy+grriUn 



-1 



(13) 



where and [/* are the radial and temporal contravari- 
ant components of the 4-velocity. The conservation law 
of baryon number can be expressed as a function of the 
proper baryon number density ub 



(«sC/^);m = g-Hg-UBU"),, 

= {nBU'),t + \{r''nBU''),r 







(14) 



The radial component of the energy-momentum conserva- 
tion law of the plasma fluid reduces to 



d_ 

dt 



dp 
dr 

-\{P + P) 



{{p + p)U'Ur) 



j,2 Qj, 



{r^ip + pp'-Ur 



dg-r 



dr 



r\2 



= . 



(15) 



The component of the energy-momentum conservation 
law of the plasma fluid equation along a flow line is 
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1 



p 



. 



(16) 



Defining the total proper internal energy density e and 
the baryonic mass density in the comoving frame of 
the plasma fluid, 

e = p-pB: pB^riBmc^, (17) 

and using the law of baryon- number conservation (p^, 
from Eq. (^6|) we have 

(eC/'');,+p(C/'');,. = . (18) 

Recalling that ^ = y(?7'^);^, where V is the comoving 
volume and r is the proper time for the plasma fluid, we 
have along each flow line 



dr 



dV 



dE 
d7 



^^^^ 



(19) 



where E = Ve is total proper internal energy of the plasma 
fluid. We represent the equation of state by the introduc- 
tion of a thermal index T(p, T) 



e 



(20) 



We now turn to the analysis of e^e^ pairs initially cre- 
ated in the dyadosphere and ionized electrons contained in 
the baryonic matter. Let n^- and ng+ be the proper num- 
ber densities of electrons and positrons associated to pairs 
created, and n^- the proper number densities of ionized 
electrons, we clearly have 



'-pair ; 



— Zrif 



(21) 



where npair is the number of e^e 
atomic number i; < Z < 1 (Z = 



" pairs and Z the average 
1 for hydrogen atom and 
Z = for general baryonic matter) . The rate equation for 
electrons and positrons gives, 



K+C/^);, 



-he 



= ai;[n^_(T)n,+ (T)-n^_n,J 



(22) 



(23) 



(24) 



where Wv is the mean of the product of the annihilation 
cross-section and the thermal velocity of the electrons and 
positrons, ne± (T) are the proper number densities of elec- 
trons and positrons associated to the pairs, given by ap- 
propriate Fermi integrals with zero chemical potential, and 



(T) is the proper number density of ionized electrons, 
given by appropriate Fermi integrals with non-zero chem- 
ical potential /ig, at an appropriate equilibrium tempera- 
ture T. These rate equations can be reduced to 



1 



crvin. 



[rie- (T)r7,e+ (T) - Tie 



- ne+ 



1 



Frac 



Equation (Ef 
(11 and li 

(r). 



n 



(T) 



ne±{T) nl_ 



(25) 
(26) 

(27) 



is just the baryon-number conservation law 
is a relationship obeyed n^i , f^e± {T) and 



The equilibrium temperature T is determined by 
the thermalization processes occurring in the expanding 
plasma fluid with a total proper energy density p, governed 
by the hydrodynamical equations (njjiajlq). We have 



P = P-t + Pe+ + Pe- + Pe- + PB, 



(28) 



where p^ is the photon energy density, ps — niBC^TiB is 
the baryonic mass density which is considered to be non 
relativistic in the range of temperature T under consid- 
eration, and /9e± is the proper energy density of electrons 
and positrons pairs given by. 



Pe± 



:Pe±{T), 



(29) 



where rigi is obtained by integration of Eq.(|25|) and 
Pe± (T) is the proper energy density of electrons(positrons) 
obtained from zero chemical potential Fermi integrals at 
the equilibrium temperature T. Whereas, p^^_ is the energy 
density of the ionized electrons, obtained by the ionization 
of baryonic matter, 

p'eAT), (30) 



pI- 



_(T)' 



where is obtained by integration of Eq.([2q) and 

Pe- {T) is the proper energy density of ionized electrons 
obtained from an appropriate Fermi integral of non-zero 
chemical potential pe at the equilibrium temperature T. 

Having intrinsically defined the equilibrium tempera- 
ture T in Eq.(p8|), we can also, analogously, evaluate the 
total pressure 

h 



P ^ P-t + Pe+ + Pe- + Pe- + PB, 

where Pj is the photon pressure, Pe± and 
r-Pe±(r), 



Pe± 



Pe 



ne± (T) 



(T) 



pUt), 



(31) 
given by, 

(32) 
(33) 



where the pressures Pe± (T) are determined by zero chemi- 
cal potential Fermi integrals, and (T) is the pressure of 
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the ionized electrons, evaluated by an appropriate Fermi 
integral of non-zero chemical potential /ig at the equilib- 
rium temperature T. In Eq.(^l|), the ions pressure ps is 
negligible by comparison with the pressures p^^f.±^f.-{T), 
since baryons and ions are expected to be non-relativistic 
in the range of temperature T under consideration. Fi- 
nally, using Eqs.( [2^j3l| ), we compute the thermal factor 
r of the equation of state . The calculation is contin- 
ued as the plasma fluid expands, cools and the e+e" pairs 
recombine, until it becomes optically thin: 



dr(n^± + ZnB)crT — 0{l), 



(34) 



where ctt = 0.665 • 10~^'*cm^ is the Thomson cross-section 
and the integration is over the radial interval of the PEM- 
pulse in the comoving frame. At this point the energy is 
virtually entirely in the form of free-streaming photons 
and the calculation is stopped. 

3. The quasi- analytical simplified model based on the 
constant-thickness approximation 

The PEM pulse expansion in the absence of bary- 
onic matter has been discussed in a previous paper 



( Ruflini et al. 1999| ) where the quasi-analytical approach 
of an expanding shell of constant thickness in the labo- 
ratory frame was adopted and validated by comparison 
with the numerical integration of the general relativistic 
hydrodynamical equations. We here generalize these re- 
sults by examining the collision of the PEM pulse with 
baryonic matter and adopting the constant-thickness ap- 
proximation both for the description of the collision and 
the further expansion of the PEM pulse by a simplified ap- 
proach to the system of equations outlined in the previous 
paragraph. 

We first recall the main results of the PEM pulse ex- 
panding in vacuum: we indicate by U{r) = Ur — const, 
the four-velocity of the slab and by 2? = r^s — r+ the 
constant width of the slab in the laboratory frame of the 
plasma fluid, the average bulk relativistic gamma-factor 7 
is, 



7 = yr+c^. 



Vr 



7 



(35) 



The evolution of the slab is governed by the total energy 
and entropy conservations, which are cast into the follow- 
ing equations as a function of the coordinate volume of 
the plasma fluid expanding from Vo to V, 



Vo 



7 = 7c 



£oVo 
eV ' 
1 dV 



(36) 
(37) 



where the proper volume V of the plasma fluid V = 7V 
and the thermal index P Eq. (pO|) , a slowly- varying func- 
tion of the state with values around 4/3, has been approx- 
imately assumed to be constant. The coordinate number 
density of e*-pairs in equilibrium is iVg± (T) = ^n^± (T) 
and the coordinate number density of e^-pairs A'^g± = 
7rte± . These equations have already been numerically in- 



tegrated (Ruffini et al. 1999) 



The baryonic matter remnant is assumed to be dis- 
tributed well outside the dyadosphere in a shell of thick- 
ness A between an inner radius rin and an outer radius 
Tout at a distance from the EMBH at which the origi- 
nal PEM pulse expanding in vacuum has not yet reached 
transparency, 



lOOrds, A = lOr, 



ds, 



A. 



(39) 



The total baryonic mass {Mb = NsmpC^) is assumed 
to be a fraction of the dyadosphere initial total energy 
(i?dya)- The total baryon-number {Nb) is then given by 



B- 



E, 



dya 



(40) 



where _B is a parameter in the range 10^^ — 10^^ and 
where mp is the proton mass. The baryon number density 
n°g is assumed to be a constant. 



Vb ' 



Pb 



(41) 



As the PEM pulse reaches the region nn < r < rout, 
it interacts with the baryonic matter which is assumed to 
be at rest. In our simplified quasi-analytic model we make 
the following assumptions to describe this interaction: 

— the PEM pulse does not change its geometry during 
the interaction; 

— the collision between the PEM pulse and the baryonic 
matter is assumed to be inelastic, 

— the baryonic matter reaches thermal equilibrium with 
the photons and pairs of the PEM pulse. 

These assumptions are valid if: (i) the total energy of 
the PEM pulse is much larger than the total mass-energy 
of baryonic matter Mb, 10"^ < E^y^/MB < 10~^, and 
(ii) the comoving number density ratio n^+^- /n^ of pairs 
and baryons at the moment of collision is extremely high 
(e.g., 10^ < n^+^-/n°g < 10^^, and (iii) the PEM pulse 
has a large value of Lorentz factor (600 < 7 < 4000). 

In the collision process between the PEM pulse and 
the baryonic matter at rout > r > rin , we impose the 
total energy and momentum conservations. We consider 
the collision process between two radii r2,ri, satisfying 
''out > r2 > ri > Hn and r2 — ri <C A. The amount of 
baryonic mass acquired by the PEM pulse is 



= -Ne± y ^ + ^- {N'^± (T) ~ N'^^), (38) AM 



Mb 47r 
Vb 3 



(42) 
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where Mb /Vb is the mean-density of baryonic matter at 
rest. The conservation of total energy leads to the estimate 
of the corresponding quantities before (with "o" ) and after 
such a collision 



(reo + Pb)7o Vo + AM = (Te + pB ^ 



AM 



'rAe)72v,(43) 



where Ae is the corresponding increase of internal energy 
due to the collision. Similarly the momentum-conservation 
gives 



(reo+p^)7o[/,°Vo = (re + pB + 



AM 



VAt)^VrV, (44) 



where radial component of the four-velocities of the PEM 
pulse \J° = \Jyi — ^ and F thermal index. We then find 



Ae = - 



(re-o + p^) 



7 



1 



7ot/°Vo 



-{VI 



- Pb -1 
AM 



AM 



(reo+p^)7oC/°Vo' 



(45) 
(46) 



These equations determine the gamma-factor 7 and the 
internal energy density e = Eq + Ae in the capture pro- 
cess of baryonic matter by the PEM pulse. After colli- 
sion (r > Tout), the further adiabatic expansion of PEM 
pulse is described by the total baryon number, energy and 
entropy conservations, i.e., the following hydrodynamical 
equati ons which generaliz e those derived in our previous 
paper ( Ruffini ct al. 199!: ) with pB ^ 0, 





V 


flB 




lo 




e 




7 


- .J- 



V7 

Vo7o ' 
r 



dt 



(iV,±) = -N, 



{reo+p%)Vo 

iT€ + pB)V ' 

V dt 



7 



(47) 
(48) 

(49) 
(50) 



In these equations (r > Tout) the comoving baryonic mass- 
and number densities are pB = Mb/V and fis = Nb/V, 
where V is the comoving volume of the PEM pulse. The 
integration is continued untill the transparency condition 
in Eq.(34) is reached. 



4. Validation of the constant-thickness approximation 

The numerical integration of the general relativistic equa- 
tions given in Sec t. have already been presented in a se - 
ries of papers (see [Wilson 1975, 1977, Wilson et al. 1996D . 
We can then proceed to compare and contrast the results 



obtained by the numerical integration of the Eqs.(14)- 
( p^ ) and the simplified quasi-analytical approach given by 
Eq.(|6|)-(||), (||)-(|4|) and (|3>-®- We select the spe- 
cific example of an EMBH of KtMq with a charge to 
mass ratio ^ = 0.1. From the equations already published 




2x10* 4x10' 
Radius/R, , .„.„ 



6x10 ' 



Fig. 1. A sequence of snapshots of coordinate baryon en- 
ergy density is shown from the one-dimensional (1-D) 
hydrodynamic calculations. This run corresponds to the 
AfBH = lO^Afo,.^ = 0.1 EMBH and B = 1.3 ■ 10''^ 



baryonic shell. The spikes are actually narrow, unresolved 
shells (see Fig. S). 




5.765 5,77 

xlO" Radius/R,,,^„^p,^, 



5.775 



Fig. 2. A detail of the last shell of the previous Figure |[ 
In this figure the thickness of the shell is resolved. 



(Prcparata et al. 1998a,t), the total energy in the dyado- 
sphere is 3.09- lO'^'^ergs. The PEM pulse is assumed to col- 
lide with a baryonic shell of thickness A = lOrds at rest at 
a radius of lOOrds- We have considered five different cases: 
(1) a shell of baryonic mass 2.23-10-''Mo (4-10^° ergs rest 
energy), from corresponding to a parameter B ~ 1.3-10"*; 
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Fig. 3. As in 

EMBH and B 



Fig. with the Mbh = IO^Mq, 
— 3.8 • 10~^ baryonic shell. 
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Fig. 5. As in Fig. |l| with the Mbh = lO^M©,^ 
EMBH and B = 3.8- 10'^ baryonic shell. 
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1000 2000 
Radius/R,,,., 



3000 



Fig. 4. As in Fig. with the Mbh = IO^Mq,^ 
EMBH and B = 1.3 • 10"^ baryonic shell. 



4000 



0.1 



(2) 6.7 • IO-^Mq (1.2 • 10 
2.2-10-3AfQ (4-1051 ergs, B 
(1.2-1052 

ergs, B ~ 



51 



ergs, B ~ 3.8 - lO"*); (3) 
- 1.3-10-3); (4) 6.7-1O-3M0 
ergs, B ~ 3.8 - lO'^); (5) 2.2 - IO'^Mq (4 - 10^2 
1.3 - IQ-^). The collision between the expanding 



slab and the baryonic matter shell is treated as an inelastic 
collision in both calculations. 

We first proceed to a qualitative analysis of the evolu- 
tion of the PEM pulse in Eqs. (14)-( p^|_f rom Figs. ^ cor- 
responding to B ~ 1.3 - lO^** and Figs7^||j| corresponding 
toB ~ 3.8 -lO""*, 1.3 -10^3^ 3.8- 10^3^ .^^Qgge i^^at the PEM 



pulse after collision with the baryonic matter shell contin- 
ues to expand as a one dimensional slab. We see, however, 
from Figs. ^ corresponding to i? ~ 1.3 - 10"^ that the ex- 




200 400 

Radius/R^y^^^^p^^^^ 

Fig. 6. As in Fig. |l| with the A/bh = lO^M©,^ = 0.1 
EMBH and B = 1.3 - 10"^ baryonic sheU. For this baryon 
shell mass we see a significant departure from the behav- 
ior seen in Figs. ^ ^, ^ |^. The baryon shell mass is sig- 
nificant enough to distort the PEM pulse shell. Thus the 
PEM pulse can no longer be reliably modeled as a shell of 
constant coordinate thickness and the comparison of Fig. 
^ can no longer be made. 
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1 10 100 1000 10 * lOS 10« 

Fig. 7. Here we see a comparison of Lorentz factor 7 for 
the one-dimensional (1-D) hydrodynamic calculations and 
slab calculations (TI/bh = 10^ A/©, ^ = 0.1 EMBH and 
B ~ 1.3 • 10"'*). The calculations show good agreement. 

pansion after collision becomes much more complex and 
the constant-thickness approximation ceases to be valid. 
From this qualitative analysis we now proceed to a quanti- 
tative analysis in Fig. ^ We compare and contrast the bulk 
gamma factor as computed from the constant-thickness 
approximation and the one from the full set of Eqs. (Q)- 
(p^). The computation refers to the case B ~ 1.3 • 10""' 
where an excellent qualitative agreement with the one- 
dimensional-slab approximation has been found. The ex- 
tremely good agreement validates the constant-thickness 
approximation. This agreement has been found up to val- 
ues of B no larger than 10~^. 

5. Considerations on the GRB structures descending 
from a constant-thickness approximation 

We now proceed to some specific prediction of GRB fea- 
tures computed by using the constant-thickness approxi- 
mation and the Eqs. (H)-®, ©-(H) and in 
the range of validation of this approximation just defined 
in the previous paragraph. As an example for clearly show- 
ing the evolution of PEM pulses colliding with baryonic 
matter, we take the following black hole mass and charge, 
as well as the mass of baryon remnant as a typical case: 



Mbh = 10-^ Mg 



i = 0.1, Mb = IQ'^E, 



dya: 



(51) 



where the total energy of dyadosphere i?dya = 3.09 • 
10^'^ergs, so the total number of e'''e~-pairs created in the 
dyadosphere (given by Eq.(|l^)) iV°+^_ = 1.9 • 10*^°, and 
baryonic mass Mb = 1.73 • 1O~^M0. This baryonic mass 
is close to the limit of validation of the slab model shown 
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Fig. 8. The Lorentz factors are given, as functions of the 
radius, in units of the dyadosphere radius, for selected 
values of shell masses given by B = (Mb(?)I E^yg, for the 
typical case Mbh = 10'' M©,^ = 0.1. The asymptotic val- 
ues 7asym = B^y^l(MB(?) = 10^, 10^, 10^ are also plotted. 



in Sect. ^. We have assumed the baryonic matter at a dis- 
tance of Tin = lOOrds, very close to the transparency condi- 
tion of the PEM pulse in vacuum (see Ruffini et al. 1999 ). 

In Fig. H we represent the Lorentz Factor of the PEM 
pulse as a function of the radius for collision with different 
amounts of baryonic matter, corresponding respectively 
to B = 10-2, B = 10-3 and B = \^-^. The diagram 
extends to values of the radial coordinate at which the 
transparency condition given by Eq.(^) is reached. Also 
represented, for each diagram, is the "asymptotic" Lorentz 
Factor: 



Tasym 



_ -Edya 



(52) 



The closer the 7 value approaches, at transparency, the 
"asymptotic" value (^2|), the smaller the intensity of the 
radiation emitted in the burst, and the larger the amount 
of kinetic energy left in the baryonic matter. This point is 
further clarified in Fig. ^, where are plotted the 7-factor at 
transparency and the "asymptotic" one as functions of the 
baryonic matter. It is interesting that, for a given EMBH, 
there is a maximum value of the 7-factor at transparency. 
After that maximum value, the energy available for the 
GRB is smaller in intensity, and at decreasing values of 
the energy, for increasing values of the baryonic mass. 

The temperature in the laboratory frame at the 
transparency point is plotted as a function of the bary- 
onic mass in Fig. |l^: it strongly decreases as the bary- 
onic mass increase. The 7T is related to the observed 



8 



Remo RufRni et al.: On the pair-electromagnetic pulse. 



1 ooooo 




1e-OOS 1e-007 1e-006 1 e-005 0.0001 0.001 0.01 0.1 

(Mb)/(Ms„|^,) 

Fig. 9. The final Lorentz-factor (the sohd hne) at the 
transparent point is plotted as function of baryonic masses 
Mb in the unit of the solar mass. The asymptotic value 
(the dashed hne) E^yi^/ {Mbc?) is also plotted. 
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Fig. 11. At the transparent point, the energy radiated 
in the burst (the solid line) and the final kinetic energy 
(the dashed line) of baryonic matter, in units of the to- 
tal energy of the dyadosphere, are plotted as functions of 
baryonic mass Mgc^ in units of solar mass. 



Fig. 10. The observed temperature Ep ^ in the lab- 
oratory frame is plotted as function of baryonic masses 
Mb<? in the unit of the solar mass. 



energy-peak of the photon- number spectrum (see e.g. 
Ruflini et al. 19991). 



We plot in Fig. |T| the energy radiated in the burst 
and the final kinetic energy of baryonic matter. We find 
that, for small values of i? (around 10~^), almost all total 
energy is radiated as GRB (see also our previous paper 
Ruffini et al. 199£), and very little energy is left as ki- 



netic energy of baryonic matter as afterglow. While for 
B ~ 10^^ roughly only 10~^ of the total initial energy 
of the dyadosphere is radiated away as GRB, and almost 
all energy is restored as the kinetic energy of the baryonic 
matter. It is also clear that for B > 10^^ the intensity 
of the Burst (see also Fig. ||) and the observed radiation 
frequency drifted to smaller and smaller values and are 
of little astrophysical interest from the point of view of 
GRBs. For such values, the energy of the dyadosphere is 
transferred practically totally to the bulk kinetic energy 
of the baryonic matter. B — 10~^ is also the limit of the 
validation of our computations based on the analytical 
slab model, described in Sect. ^. For values of baryonic 
matter in between these two extremal cases {B ~ 10~^ 
and B = 10~^) the analytical slab model covers the whole 
range of the observed properties of Gamma Ray Bursts. 

We turn now to the thermodynamic parameters rel- 
evant in the evolution of the PEM pulse. In Fig. ^ the 
temperature of PEM pulse, both in the comoving and in 
the laboratory frame, are given as a function of the radius 
for a typical case (Mbh = lO^Mg, ^ = 0.1 and B = 10~^). 

In Fig. |l^, we plot the total energy of the non bary- 
onic components of the PEM pulse, which includes both 
thermal and kinetic energy, and the kinetic energy of the 
baryonic matter as functions of the radius, for the typical 
case Mbh = lO^M©, ^ 0.1 and B ^ 10"^. The total en- 
ergy of the non baryonic components of the PEM pulse is 
equal to E'dya before the collision (see Ruffini et al. 1999 ) 
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Fig. 12. The temperature in the comoving frame r'(MeV) 
(the solid line) and in the laboratory frame 7T' (the 
dashed line) are plotted as functions of the radius in the 
unit of the dyadosphere radius r^s for the typical case 
Mbh = lO^Mo,^ = 0.1 and B = 10'^ (see Sect. |in the 
text) . 
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Fig. 13. The total energy of the non baryonic components 
of the PEM pulse (the solid line) and the total kinetic 
energy of baryonic matter (the dashed line) in the unit of 
the total energy are plotted as functions of the radius in 
the unit of the dyadosphere radius r^g for the typical case 
Mbh = 1Q^Mq,1^ 0.1 and B = 10"^^ 



and drops after the collision. While the kinetic energy of 
baryonic matter 



Ek - PbV{^ - 1) 



(53) 



increases as function of radius for (r > ri„). The sum 
of both them is equal to the total energy i?dya — 3.09 • 
lO^^'ergs during the evolution of the PEM pulse. The value 
of the total energy of the non baryonic components of the 
PEM pulse at the transparency point, the ending point of 
the curve in Fig. |l^, is the energy released in the burst. We 
have discussed this energy as function of baryonic masses 
in Fig. |ll|. 

Before and after the collision, the condition of entropy 
conservation applies, and we have: 



■Sbcforc =^ {P'f+ Pe± +P7 +P(e±)) 



''after ^ 



+ Pe± + Pe- + 



(54) 



(55) 



where e_B is the thermal energy of baryonic matter, and 
we can neglect, in Eq. (|5|), the pressure of the baryonic 
matter. A sudden increase of the entropy occurs during the 
collision both for the addition of the baryonic matter, and 
for the thermal reheating due to the inelastic collapse of 
the PEM pulse with the baryonic matter at rest. From the 
energy and momentum conservations, we obtain for the 



values of the 7 and the temperature during the collision: 
the proper internal energy, i?int, of the plasma increases 
as 



dEi,,t - (7 - 



B 



(56) 



and the slab is decelerated, in terms of Lorentz factor 7, 
as 



(i7 — 



7 



1 



Mbc2 + Ei, 



-dMB , 



(57) 



as baryon mass, cIMb, is incrementally gained. 

Before the collision the PEM pulse expands keep- 
ing its temperature in the laboratory frame constant, 
while its temperature in the comoving frame falls (see 
[Ruffini et al. 1999| ). During the collision, a heating of the 
plasma due to the energy and momentum conservation oc- 
curs (see also Fig. |l^, where reheating process leads to an 
increment of the number density of e"'"e~ pairs). As the 
system expands further, both the comoving temperature 
and the temperature in the laboratory frame decreases, 
since the total energy of the e+e" pairs and the photons 
before the collision is constant and equal to E^ya, while 
after the collision 



E, 



■dya 



Ep 



E,. 



E, 



photons 



(58) 



where i?Baryons is only the thermal energy of the baryons. 

It is also interesting to monitor the change of tem- 
perature in the comoving frame before the collision and 
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Fig. 14. The number densities Ue+e- (T) (the solid line) 
computed by the Fermi integral and n^+^- (the dashed 
line) computed by the rate equation ( |38| ) are plotted as 
functions of the radius. T' <C irieC^, two curves strongly 
divergent due to e^e~-pairs frozen out of the thermal equi- 



librium. The peak at r = 
reheating in the collision. 



is due to the temperature 



after the collision (see also Fig. [T^ ). Before the colli- 
sion, due to the entropy conservation, in the process of 
e+e~ annihilation the factor r^-^y (where the subscript 
"o" means "initial value") increases to a value near to ^ 



(see Weinberg 1972 ; Rufhni ct al. 1999 ) since the collision 
occurs at r = lOOrds, near the condition of transparency. 
The number of e~^e~ pairs has now been reduced drasti- 
cally. The further jump in the value of the ratio ^^tt- is 
principally due to the energy and momentum conservation 
during the inelastic collision. After the collision, there is a 
small reheating due to the annihilation of the e^e^ pairs 
created in the collision (sec Fig. |l^, which cannot be seen 
on the scale of our plot. Then, the ratio remains constant. 

6. Conclusions 

By the direct comparison with the numerical integration 
of the complete relativistic hydrodynamical equations, the 
use of the constant-thickness approximation has been vali- 
dated for values of the parameter B < 10~^. For B > 10~^ 
the amount of energy released at transparency in the burst 
decreases, and its energy drifts toward low energy values, 
which are of little interest for the astrophysics of GRB. 
We conclude that the constant-thickness approximation 
is valid for all astrophysically relevant situations for the 
analysis of the GRB at transparency. Based on this vali- 
dation we have studied the evolution of a PEM pulse in 



Fig. 15. We plot the ratio {T^V)/ {T^Vo), T being the 
temperature of the PEM pulse, and V the volume in the 
comoving frame, as function of the radius in the unit of 
the dyadosphere radius r^s for the typical case A/bh — 
IO^Mq,^ = 0.1 and B = 10"^ (see Sect. | in the text). 
After the collision, there is a small reheating up to the 
annihilation of the e'^e~ pairs recreated in the collision 
(see Fig. 14), which cannot be seen on the scale of our 
plot. 



the presence of selected amounts of baryonic matter. We 
have studied for a typical case of Mbh = lO^M©,^ = 0.1 
EMBH and B = 10"^ baryonic shell, the thermal, the 
bulk Lorentz factor evolution of the PEM pulse as well as 
the kinetic energy left over in the baryonic matter after 
the decoupling of matter and radiation and the emission 
of the GRB. Additional results corresponding to a larger 
range of masses and charges of the EMBH and the cor- 
relations between the peak energy and the duration of 
GRBs to be expected in our model will be considered in 
a future paper, together with the results of analyzing the 
interaction of the kinetic energy left over in the baryonic 
matter, after the decoupling of matter and radiation, with 
the interstellar medium. 
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